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Abstract 

Current induced switching in Co/Cu/Co trilayers is described in terms of ab-initio determined 
magnetic twisting energies and corresponding sheet resistances. In viewing the twisting energy as 
an energy flux the characteristic time thereof is evaluated by means of the Landau-Lifshitz-Gilbert 
equation using ab-initio parameters. The obtained switching times are in very good agreement with 
available experimental data. In terms of the calculated currents, scalar quantities since a classical 
Ohm's law is applied, critical currents needed to switch magnetic configurations from parallel to 
antiparallel and vice versa can unambiguously be defined. It is found that the magnetoresistance 
viewed as a function of the current is essentially determined by the twisting energy as a function 
of the relative angle between the orientations of the magnetization in the magnetic slabs, which 
in turn can also explain in particular cases the fact that after having switched off the current the 
system remains in the switched magnetic configuration. For all ab-initio type calculations the fully 
relativistic Screened Korringa-Kohn-Rostoker method and the corresponding Kubo-Greenwood 
equation in the context of density functional theory are applied. 

PACS numbers: 75.30.Gw, 75.70.Ak, 75.70.Cn 
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I. INTRODUCTION 



Reversal of the orientation of the magnetization without applying an external field seems 
to be of considerable interest for magnetic switching of micro-devices and caused extensive 
experimental and theoretical studies of the effect of currents on magnetic nanostructures. 
The experimental fact a ll2l3l4l5l6l7l8l9ll ° l11112113 are still quite confusing, or, to put in the words 
of a recent short review article 1 ^, "many observed phenomena can be described qualitatively 
... by a simple semi-classical spin-torque model. However, evidence of complications from 
several experiments suggests that a full understanding of all observations is not yet achieved" . 
The by now generally accepted experimental facts are the following ones: (1) if the current 
of a given sign favors the parallel (P) magnetic configuration, the current of the opposite sign 
favors the antiparallel (AP) configuration, (2) the current needed to switch the magnetic 
configuration in nanostructured magnetic multilayer systems is of the order of 2 - 5 mA 
in samples with a volume in the range of 40 - 800 nm 3 . These two facts led inter alia to 
a schematic effective two-level energy diagram for switching in which the critical current 
corresponds to the energy needed to overcome the potential barrier between the parallel 
and the antiparallel magnetic configuration. Since experimentally also telegraph noise is 
observed, which in turn seems to correspond to an oscillation between these two states, 
this schematic picture proved to be quite useful. If by means of a sufficiently high current 
the system is driven from one configuration to the other one, it also can happen that after 
turning off the current the system remains in the switched configuration, i.e., the system 
does not return to the ground state. In the two-level energy diagram this would correspond 
to the case that the two schematic minima are separated by a high enough potential barrier 
and are of about the same energy. 

Most theoretical investigation o 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 were concerned about finding ex- 
pressions for the interaction between the applied current and the orientation of the magnetic 
moments. Almost all theoretical considerations and models used the concept of spin currents 
and had to use phenomenological parameters to relate the respective approach to the exper- 
imental evidence. Quite clearly in most cases the main idea was to describe the cause for 
current induced switching and deal afterwards with the subsequent effect, namely a change 
in the magnetoresistance. Therefore the effect - creating excited states - was interpreted 
in various ways by invoking spin waves, all kinds of spin-polarization effects, etc. It is be- 
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yond the scope of the present paper to summarize the various theoretical approaches used 
up-to-now. 

In here a completely different approach is pursued: the main idea is to calculate fully 
relativistically the twisting (exchange interaction) energy of a system as it goes from a 
parallel to an antiparallel configuration, or opposite. This is a continuous function of the 
relative angle between the orientations of the magnetization in the magnetic parts of a spin 
valve system. In keeping one orientation fixed and rotating the other one by an angle 
around an axis perpendicular to the fixed orientation one thus can switch continuously from 
say the parallel magnetic configuration to the antiparallel configuration. For each given 
rotation angle simultaneously the corresponding sheet resistance (resistance divided by 
the unit area) is calculated fully relativistically, which then is also a continuous function of 
the rotation angle. It should be noted that by using a fully relativistic approach the spin 
no longer is an observable, i.e., at a given angle there is just one sheet resistance. In 
adopting this approach (1) the excitation energy is related to the rotation angle, and (2) for 
the same angle a physical observable, namely the sheet resistance is evaluated. Therefore at 
a given the effect of the physical phenomenon is described, which then can be related to 
the cause, namely the turning on of a current. It will be shown later on that the twisting 
(exchange interaction) energy is the ab-initio analogon of the above mentioned two-level 
energy diagram for switching. Furthermore, by means of relating the twisting energy and 
the corresponding sheet resistance to the current not only a critical current can be defined 
unambiguously, but also the complexity of the switching process becomes evident. Quite 
clearly in this picture no dynamic effects can be calculated, although very good reasons for 
the occurrence of the telegraph noise can be given. The quantum mechanical tools applied 
are the fully relativistic screened Korringa-Kohn-Rostoker method^ and the fully relativistic 
Kubo-Greenwood equation^! in the context of the local density functional approximation. 
All further reasoning is based on the Landau-Lifshitz-Gilbert equation, in terms of which 
switching times can be evaluated using ab-initio parameters. The introduced approach is 
applied to Co/Cu/Co type spin valves and in fact will show quite a few of the experimentally 
observed features mentioned earlier. 
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II. TWISTING ENERGIES AND SHEET RESISTANCES 



Consider a typical trilayer system of the type FM/NM„/FM consisting of two semi- 
infinite magnetic leads (FM) and a so-called non-magnetic spacer (NM) such as for example 
Co(100)/Cu„/Co(100) or equivalently Co(100)/Cu n /Co m /Vac, where n denotes the number 
of spacer layers and m is a sufficiently large number of layers of the magnetic metal. Sup- 
pose now that ~uq denotes a particular unit vector (reference orientation, either parallel or 
perpendicular to the surface normal) characterizing the orientation of the magnetization in 
a particular atomic layer containing one atom per unit cell. If 

~riB = ~nB' = ~ni = ~no , i = l,...,n, (1) 

such a configuration is usually referred to as a parallel configuration, whereas for 

~nB = ~ni = ~no, i = 1, • • • , (n/2) ; ~n B , = ~n i = -~n , i = (n/2) + 1, . . . , n , (2) 

frequently the term "symmetric" antiparallel configuration is used. If ~n b-, ~n b 1 (the ori- 
entations of the magnetization in the semi-infinite leads) and the ~n \ are each rotated by 
individual angles around an axis perpendicular to ~n* this situation refers to a general non- 
collinear magnetic configuration in two-dimensional translational invariant systems. As for 
reasonably large n the interior of the NM part is completely non-magnetic in the following 
specific non-collinear configurations of the type 

~nB = ~ni = ~no , i = 1, . . . , (n/2) ; ~n b> = ~n i = ~n' , i = (n/2) + 1, . . . , n , (3) 

will be considered, where n' is a unit vector rotated by an angle 9 with respect to Iiq. 
It is obvious that for these magnetic configurations it is sufficient to specify the rotation 
angle 6. Expressed in simple terms this means that in the right half of the trilayer system 
the orientation of the magnetization is rotated uniformly by an angle with respect to the 
orientation of the magnetization in the left half. 

Since a current perpendicular to the planes of atoms has to be described in the present 
paper the reference orientation 7t is chosen to be parallel to the surface normal (z axis); 
the rotations are performed around the y axis. 
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A. Twisting energies 



The energy difference between the two possible collinear states, namely the parallel and 
the "symmetric" antiparallel magnetic configuration, is usually termed interlayer exchange 
coupling energy. In using the magnetic force theorem 2 ^ the total energies of these states are 
replaced by the corresponding grand potentials (at zero temperature), i.e., by the so-called 
band energy difference AE®(it), 

A£ 6 °(tt) = E b (n) - E b (0) = E b (AP) - E b {P) , (4) 

j-E F 

E b (Q)= n(0; E)(E — E F )dE , (5) 

Je 

where n(0; E) is the density of states for a particular configuration, E Q the valence band 
bottom and Ep the Fermi energy. According to Eq. (J1J the below convention applies 

n > P : ground state 

AE°M = \ g . (6) 

I < AP : ground state 

In a similar manner the "twisting energy" is defined by the following difference, 

AE b (Q) = E b (Q) - E b (0) , < < vr , (7) 
and - as is well-known - can be expanded in a power series in cos O 

AE b (&) = a(l - cos(e)) + 6cos 2 (6) + ccos 3 (9) + . . . , (8) 
such that in all orders 

a = AE b (n/2) . (9) 
In first order AE b (Q) is then approximated by 

AE b (Q) ~ AE^iQ) = a(l - cos(6)) , (10) 

in second order by 

AE b (Q) ~ AEf } (6) = a(l - cos(6)) + 6cos 2 (6) , b = AE b (n) - 2AE b (n/2) , (11) 
in third order by 

AE b {Q) ~ A^ 3 \Q) = a(l - cos(6)) + 6cos 2 (6) + ccos 3 (6) , (12) 



b = - AE b (ir) - 2AE b (n/2) + 8AE b (2n/3) , c = 8AE b (2n/3) - 2AE b (n) , (13) 

etc., where AE b (7r), AE b (ir/2), . . . , refer to the actually calculated values. 

It should be noted that most frequently AE b (&) ~ AEj ) 1 \&) is assumed, an approxima- 
tion, which - as will be shown later on - not necessarily is granted. Clearly enough in using 
only A.E^(0) the coefficient a is simply half of the interlayer exchange coupling energy, 
a = A-Ej (7r)/2. In principle by calculating AE b (Q) for a few selected values of O, in terms 
of Eq. (jHJ) a reasonably good approximation to AE b (Q) for 6 varying continuously between 
and 7i can be obtained. 



B. Sheet resistance and magnetoresistance 

As is well-known in CPP (current perpendicular to the planes of atoms geometry) the 
magnetoresistance can be defined via the sheet resistances for the respective collinear mag- 
netic configurations P and AP, 

Mii W = ^M ^(AP)-r(P) 

v ' r(vr) r(vr) r(AP) v ' 

since the resistance R(Q) is defined as 

R(S) = r(e)/A , (15) 

where Aq is the unit area. In a similar manner for the present non-collinear configurations, 
the difference in sheet resistances is given by 

Ar(6) = r(0) - r(0) , (16) 

and the corresponding magnetoresistance by 

Mfi(O) = *§> . (17) 

The difference in sheet resistances can again be expanded in a power series in cos 

Ar(0) =«(l-cos(6)) + /3cos 2 (e)+7cos 3 (e)+ ... , (18) 

with a = r(7r/2) — r(0). It will be shown later on that in most of the cases investigated in 
here 

Ar(0) ~ Ar^(0) = a{l - cos(6)) , (19) 



i.e., 



r(G) = r(0) + a(l - cos(G)) 



(20) 



C. Magnetic Joule's heat generated by a current 

For I — 0, G takes on its equilibrium value O eq . As the current / is turned on, the 
relative orientation of the two magnetic layers changes to G. Evidently, the work done to 
accomplish this rotation is AE(Q) = E(Q) — E(Q cq ). Suppose that this energy difference 
is equal to the energy lost by the current in the form of a "magnetic" contribution to the 
Joule's heat Q, 

Q = R(Q)I 2 . 

Thus for a fixed current /, 

AE(S) = ri?(G)/ 2 , 

where r is the time required to accomplish the rotation. This equation can be solved for the 
function Q (J), whose inverse is given by 

/(G) = iv^v^epA^e) , (21) 

where 

AE(Q) = E(Q) - E(0) + min [AE(Q)\ , (22) 

i.e, where AE(Q) is a positive definite excitation energy. 

Since Q, AE(Q) and therefore /(Q) are scalar positive definite quantities, the above 
construction is independent of the direction of the current flow. Nevertheless, in the follow- 
ing, the concept of twisting energies AE(Q) and the corresponding magnetic Joule's heat 
generated during a time interval r shall be used to explore the physics of current induced 
switching. In short, evidence will be provided that the origin of the work done against the 
exchange forces acting between the two magnetic layers is the magnetic contribution to the 
energy dissipation from the current. 

It might seem that by using an energy flux relation the problem of evaluating the current 
I was only shifted to yet another unknown quantity, namely to the characteristic time r, 
whose theoretical description and evaluation therefore has to be the subject of the next few 
sections. Furthermore, it has to be pointed out that any comparison with experimental data 
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has to take into account also the actual area Aq present in a given experiment. However, 
before going ahead to discuss these two quantities, the computational details of the ab-initio 
related parts of this paper shall be given. 

D. Computational details 

The effective scattering potentials and exchange fields of spin valve systems of the type fcc- 
Co(100)/Coi2/Cu n /Co m /Co(100), 12 < n < 36, m > 11 were determined selfconsistently 
using the fully relativistic spin-polarized Screened Korringa-Kohn-Rostoker method,— where 
at least m layers of Co served as "buffer" to the semi- infinite leads. It should be noted 
that because of the special features of the applied screened structure constants^ the total 
number of atomic layers between the two semi-infinite systems has to be a multiple of three. 
For this reason the thickness of the right buffer had to be kept variable. In all cases the 
local density approximation of Vosko et al.-^ and, in order to obtain selfconsistency, a total 
of 45 k\\ points in the irreducible part of the surface Brillouin zone (IBZ) was applied. All 
selfconsistent calculations were performed with the orientation of the magnetization pointing 
uniformly perpendicular to the planes of atoms (reference configuration). 

In using the magnetic force theorem the twisting energies with respect to this reference 
configuration were then evaluated for each n using the symmetric arrangement, i.e., for 
the left half of the system (fcc-Co(100)/Coi2/Cu n /2) the orientation of the magnetization 
remained unchanged whereas the right half was rotated by a particular (uniform) angle G. 
For this kind of calculation a total of 960 kn points in the IBZ was used, a set-up, which yields 
very reliable results. For further computational details concerning the Screened Korringa- 
Kohn-Rostoker method and the evaluation of band energies, see the review article in Ref. 

y 

The sheet resistances for a given rotation angle G were first evaluated by means of the fully 
relativistic version of the Kubo-Greenwood equation at Ep + i5, 5 > 0, and then numerically 
continued to the real energy axis. In this part of the calculations for the occurring Brillouin 
zone integrals a total of 1830 k» points was used. For a detailed discussion of this approach, 
see the review article in Ref. ^71 
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III. THE LANDAU LIFSHITZ GILBERT EQUATION 



From the polar form of the Landau-Lifshitz-Gilbert (LLG) equation,— 
dM 7g 



dt 1 + a G l 



-M x H cS + ^-M x (m x H eS 



( 7G ,a G >0) , (23) 



where 7g is the Gilbert gyromagnetic ratio (precession constant), «g is the dimensionless 
Gilbert damping parameter, M is the magnetization, Mo = M , and H eS is the local 
and time-dependent effective field, one immediately observes that (1) M precesses almost 
purely, if damping is low («q ~~ > O), 31 (2) almost no precession, but slow switching occurs, 
when the damping is high (oq — > oo) 32 , and (3) the fastest switching refers to oq = l. 33 
Rewriting the LLG equation in terms of an experimental damping parameter G,— 

1 dM .-s - s G ( - dM 



-M x # + [M x -1- , (24) 



7g dt " 7g-^o V ^ 

the dimensionless Gilbert parameter «g is given by 3 * 5 - 



«G = — %-r , (25) 

and the Gilbert gyromagnetic ratio 7g b y 17131136137 

9PB 9\e\ , * 

7G = ~ =2^' (26) 

where e refers to the elementary charge, m e to the mass of an electron, /ib to the Bohr 
magneton, and g is the (electronic) Lande-factor. Experimentally , 35 ! 38 it has been shown 
that in multilayer systems G in Eq. (|25j) varies linearly with 1 / d, where d is the film 
thickness, see also Tab. [I] 

By definition the magnetization M refers to the volume averaged total magnetic moment. 
Assuming, however, that in a layered system the layer-resolved magnetic moments Mj, 
where % denotes atomic layers, are coherently precessing,^ it is sufficient to describe the 
magnetization dynamics of the layered system in terms of the motion of either the layer 
averaged magnetic moment M, 

f --7Mx#- + Q |x(Mx^), (27) 

1 - 

i=l 



where N denotes the number of magnetic layers, e.g. N = m + n / 2, or in terms of the 
magnetization direction n ; 30 i 40 

''" x H cS + alt x (ltxH cS ) , n=^. (28) 



dt 1 V / ' M f 



'0 



It should be noted that the equivalence of these two equations, namely Eqs. (|2*7j) and (J2*%|) . 



relies on the conservation of M. , which in turn implies that M 



Mq and | n \ = 1. 



A. Internal effective field 

The local effective field H eG that enters Eqs. ()27|) and (|28|) can directly be derived from 
the Helmholtz free energy density by taking its variational derivative with respect to the 
magnetization , 33 i 35 i 36 i 41 



^=-Sr- v ^' with T=F v (29) 



where V is the total characteristic volume of the system and the free energy F includes the 
exchange energy, the crystalline anisotropy energy, external magnetic fields, etc. , 1 either 
in a parameter-free manner or by using different types of model Hamiltonians. Since for 
layered systems as considered in here, 

d 



Eq. (|29|) can be written as 



M=x,y,z p 



H cS = ^ = -V^F , (30) 
DM M 



where the Helmholtz free energy F, 

1 - 

i=i 

refers to the reference volume ^ and is again the number of magnetic layers considered. 

According to Eq. ()30|) the internal effective field, H E , arises from the contribution of the 
total energy E b to the free energy, 

H E = — — = —VjtfEb . 
dM 
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Assuming, e.g., in terms of Eq.(JHJ) that the derivatives 

}fc TP. ( tut . \ 3 



8MtdMy 2 dMt . l 

in the below Taylor series expansion of the total energy, 

E b (M + M) ~ E b (M ) + J] 1 (M ■ V Al ) " E b (M ) , (31) 

k=l 

are available up to a certain order p, where M is the initial reference moment, the Cartesian 
components of the internal effective field can directly be given. 

In particular, if the change in the moment M is constrained to the Oyz-plane, 

M = M y e y + M z e z , (M x = 0) , 

then 

(32) 



^e_ y-y- d k E(M ) 



k=i g=o um y olvIz 



y z — * i y z ^ 

re y + — ^-ttt —e z 



(k-q-l)\q\ y (k-q)\(q-V 

with Mq = Mqc z being the initial, ground state moment. Provided that the magnitude of 
the moment is preserved, 

Ml = M 2 + Ml = M 2 , (33) 

by keeping in Eq. (|3Tj) only terms up p = 3, one gets 

/ -\ - - M M 2 M 3 

AE b \M) = E b (M + M) — E b (M ) ~ a - + 6 ^ + , (34) 

where the coefficients a, 6 and c are defined in Eqs. © and (|13j) . Therefore the energy 
torque rotating the moment is given by 

M x H E = e x M y H E = -n y (-a + 26 n z + 3cn 2 z ) e x , (35) 

whereas 

M x (m x = -?iy (-a + 26 n z + 3cnf) (M z e y - M y e z ) . (36) 
B. The characteristic time of switching 

Inserting the internal effective field H E into the LLG equation (in the absence of precession 
around the z-axis), 



dM M t -* -* p \ 7g 
-«irr x \MxH E ), a = a G 



dt M V " ' / 1 + a. 



2 ' 
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then leads to 



Mo— ^ = 
at 



Mi 



dt 
dn z 



-an y n z 



(—a + 2bn z + 3cn^) 



an' 



(—a + 2bn z + 3cn 2 ) 



(37) 
(38) 
(39) 



Since according to Eq. (j3Hj) . n 2 + n 2 =1, Eq. (J3~9"j) reduces to 



M 



a 



(l-nl) (-a + 2bn z + 3cn 2 z ) 



(40) 



Assuming that n z ^ ±1 or (—6 ± vP + 3ac) / 3c and b 2 + 3ac > 0, c ^ 0, Eq. (jlUj) can 
directly be integrated and leads to the time r = tf — U needed to change n z from n z = n z (tj) 
to n{ = n z (t f ) 



a 

M 



1 



-T 



2 [(3c - a) - 2b] 
b 



In 



7? 



ni + l 



2 [(3c - a) + 2b] 



In 



7? 



£-1 



nl - 1 



+ 



(3c - a) 2 - 46 2 
1 



In 



3c (n{) + 2bn f z -a 



3c (nl) + 26 Tij - a 



a (3c - a) + 26 2 



(3c -a) -Ah 2 2VW+3ac 



x In 



(b + 3cn() - Vb 2 + 3ac (6 + 3cn») + y/b 2 + 3ac 



(6 + 3ct7* ) - V6 2 + 3ac ^ + 3cn / j + ^yi + 3ac 



(41) 



Thus the time r (ag) needed to change n z from 77* to 77^ as a function of the Gilbert damping 
parameter oq > 0, when using a = q;g7g/ (1 + c^g) an d Eq. PBjl . can be written as 



T(aG) = c ^i±^ = ft ^ c i±^ 



(42) 



7g «G # «G 

where C denotes the rhs of Eq. (j41|) and DJIq is the magnetic moment M in units of Bohr 
magnetons /xb- 

According to Eq. (|42j) it follows that the minimal time of changing n z from 77* to n( is 
given by 



f 1\ OK m ° n A t \ i + ^G^iin 

r m i n = r (a G = 1) = 2/7 C and r (a G ) = — 

9 a G 2 



(43) 



12 



IV. RESULTS 



A. Magnetization and switching times 

The selfconsistently obtained constant magnitude of the magnetization for fee 
Co/Cu/Co(100) and in the case of a thin Co slab of M f cc Co = 1.418 x 10 6 A m" 1 is 



in very good agreement with the available experimental data for Co bulk, see Ref. 



43 or 



441 . By using the experimental Lande g-factor for fee Co, namely g icc Co = 2.146 ± 0.02^ 



7 fcc Co = 18.87213449 x 10 10 m A -1 s -1 , which in turn yields the below ratio 

]\/ficc Co fYTJ-fcc Co 

f ° r = h ° = 0.782531193 x 10" 34 Js = 4.884173446 x 10~ 16 eVs . 

-yg c Co gfee Co 

Since according to Tab. [n]the quantities on the rhs of Eq. (|41|) are of the order of (meV) 1 = 
10 3 (eV) -1 , this implies that the time needed to change the moment direction from n* to 
n[ is of the order of 10~ 13 s, namely femtoseconds. From Tab. IIHI one immediately can 
see that the theoretically obtained values of r m i n are within the range of values known from 
micromagnetic simulations for a polycrystalline thin Co film, which showed that the reversal 
time ranges from 0.05 ns for «g = 1 to 0.2 ns for (Xq — 0.1 . 43 It should be noted that because 
the sign of r min is uniquely determined by the sign of the initial and final values for n z , in 
Tab. Illll onlv those values for n\ and n( are listed, which yield r min > 0. The sign of the such 
determined n( confirms therefore independently the ground state configuration predicted 
by the magnetic force theorem, see Eq. (J7J). 

As can be seen from the corresponding column in Tab. IIHI the switching time is largest 
for the Gilbert damping parameter of Co bulk. By scaling «g to the thickness of the 
thin (rotated) Co slab used in the present calculations according to values found for 
Coat/Cu/Co(100), 44 the magnetization reversal time r (a G ) in Eq. (|43[) changes only very 
moderately in comparison with r m i n . 

B. The importance of cross sections (unit areas) 

Going now back to Eq. (j2"Tj) . rewritten below by indicating the appropriate units, 



7(e)i -™ 7i437 -^V^Sb mA ' (44) 
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it is obvious that for any kind of comparison to experiment not only r has to be evaluated, 
but also that Aq, the surface perpendicular to the z axis through which the current / (0) 
flows, has to be taken into account. Usually the cross section of nanopillars is given in 
nanometer (nm), i.e., is of the order of 

(^o)si = (Ao) nm i x 10- 18 , 
which combined with the switching time (in nanoseconds), 

(r)si = (r> M x 10- 9 , 

yields the following factor that multiplies the square root of the (quantum mechanically 
derived) quotient of twisting energy and sheet resistance in Eq. 



'<4>>si_. /<A>> nm2 xlO- 18 __ l(A ) nm2 



x 10" 9 



(t)si V H,xW-" V H 
In using, e.g., (A ) nm2 = 120000^, other values of (A ) nm2 are listed in Tab. IIV[ and 
(Tmin) ns = 0.01, this results into a value for a/ (^o)si / ( r min) SI of about 0.11. 



C. Twisting energies and currents 

In all figures showing twisting energies and sheet resistances, etc., the actually calculated 
values are displayed, solid lines only serve as guidance to the eye. For illustrative purposes 
also the first order approximation to the twisting energy is depicted in these figures as a 
dashed line. As it is not possible to show all results obtained these figures concentrate 
on systems with the spacer thickness varying between about 35 - 50 A. This still results 
in a considerable number of figures, which, however, seems to be necessary considering 
that in experimental studies mostly nanopillars are used, i.e., most likely an average over 
thicknesses is recorded, and also in order to illustrate the complexity of the effects to be 
seen. Furthermore, in all figures for the current / (6) the factor */ (A ) SI / (T m i n ) SI in Eq. 
(JUJ) is replaced by unity. 

In the investigated range of spacer thicknesses the number of cases in which the twisting 
energy is proportional to (1— cos(G)), see the corresponding figures for n > 31, is surprisingly 
small, whereas in all cases the sheet resistance - more or less - is of this shape. This in 
turn implies that all special features to be seen for are mostly related to the functional 
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form of the twisting energy. Taking for example n = 20, 1(0) remains about constant for 
> 90°, a value which refers also to the critical current that has to be applied to drive the 
system from parallel to antiparallel. However, one also can see from this figure that AE(,(Q) 
has a maximum at about 140°: the system has to overcome a small barrier to return to 
the ground state (parallel configuration). From the entry showing the magnetoresistance 
versus current, it is evident that at the critical current the magnetoresistance jumps by 
about 20%. For n = 21 the situation is even more dramatic, since AEf,(0) has quite a large 
maximum at 90°, the AP configuration being only slightly less energetically favored than 
the P configuration. In this particular system the meaning of the critical current is quite 
obvious: it simply is the maximum in the /(0) versus curve. The same situation, even 
more impressive, pertains for n = 25, since now the parallel and antiparallel configuration 
are virtually degenerated in energy, separated, however, by quite a barrier. The figures for 
n = 20, 21, 23, 24, 25 and 27 are perfect ab-initio analoga for the schematic effective two-level 
energy diagram mentioned in the introduction: the energy displayed in this ad hoc scheme 
is nothing but the twisting energy, the schematic abscises being the relative angle between 
the two orientations of the magnetization. 

The system with 26 Cu spacer layers is in particular interesting, since a non-collinear 
configuration is the ground state. In this case a tiny current (about 0.05 mA in this figure) 
produces a magnetoresistance that can be either zero, 1.5 or about 6 %. It should be noted 
that the energy barrier between the ground state and the parallel configuration is minute: 
the system can almost freely oscillate between magnetic configurations for values of < 60°. 
Another interesting case seems to be for 30 Cu spacer layers, which shows a strong deviation 
from A£ 6 (1) (0) at about = 90°, not enough, however, to cause an additional minimum 
between the AP (ground state) and the P configuration. For < < 100° the current 1(0) 
varies almost linear, changes slope, and varies again almost linear for > 100°. 

For n > 31 no more interesting effects are observed: the twisting energy can be described 
very well in terms of AE^'(0); in order to switch from parallel to antiparallel or vice versa 
a current of about 0.35 - 0.6 mA is needed. It should be recalled that all values of 1(0) 
quoted in this section refer to *J (A)) si / ( r min)si = 1 in Eq. (jH)) . 

In the last figure finally the rotation of the magnetization around the z-axis (precession) 
is shown for n = 25. As can be seen the precessional changes in the twisting energy are very 
small indeed. This figure justifies a posteri the approach taken to evaluate and discuss the 
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switching time r min . 
V. DISCUSSION 

In viewing now all the various cases discussed above the following observations can be 
added: (1) if the slope of the magnetoresistance with respect to the current is uniformly 
positive (negative), the parallel (antiparallel) magnetic configuration is favoured (see, n = 
31,33 versus n = 32), (2) if it becomes approximately infinite at a certain current then a 
jump in the magnetoresistance occurs (e.g., n = 20,24,28), and (3) if this slope changes 
sign, a more complicated behavior pertains (e.g., n = 21,22,25,26). In the latter case the 
system either remains in the switched configuration (n = 21, 25) or because of a non-collinear 
ground state oscillations in the magnetoresistance between zero and a few percent can occur 
when a very small current is applied (n = 26). The so-called telegraph noise seems to refer 
to the jumping between such minima in the twisting energy, the jumping rates obviously 
being connected with the barrier between these minima. The current needed to switch a 
configuration from parallel to antiparallel (or vice versa) refers to the largest value of /(©). 

The present results suggest that the efficiency of current- driven switching can considerably 
be optimized by varying the spacer thickness: theoretically in using a spacer thickness of 
about 43 A (25 layers of Cu) perfect switching can be achieved such that the system remains 
in the switched state after the current is turned off. Further theoretical investigations using 
the approach presented in here can include interdiffusion effects at interfaces or refer to 
different kinds of magnetic slabs (leads) such as for example permalloy, a system, which 
because of anisotropy effects in the magnetoresistance perhaps is even more complicated 
than the present Co/Cu/Co trilayer. 

Altogether correlating the twisting energy and the corresponding resistance with the 
current yields a very consistent view of the complexity found in current-driven experiments. 
Clearly enough this correlation suffers from the fact that up-to-now no quantum mechanical 
description for the Gilbert damping factor was found and that a linear response theory 
(Kubo-Greenwood equation) is used to evaluate the electric transport properties, i.e., that 
the current had to be formulated as a scalar quantity. 

Finally, it has to be remarked that the experimentally observed critical switching cur- 
rents are by a factor of about 10 - 100 larger than the ones obtained in here. It has to 
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be remembered, however, that in here an ideal Co/Cu/Co trilayer was assumed while most 
experiments are based on rather complicated nanostructures such as for example nanopil- 
lars and therefore - although a consistent approach to was introduced in order to evaluate 
switching times in terms of ab-initio parameters - also the question of a comparable cross 
section (unit area) is of quite some importance. 
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TABLE I: Experimental damping parameter G for different systems. 



material 


type of system 


G (10 8 s- 1 ) 


Fe 


bulk 


0.52i, 0.58^ 






0.59 ± 0.06; 0.572 ± 0.04 45 






0.8*1 0.7 ±0. 06^ 


Fe 


single film 


I 5 35.38 






1.3±0.li£ 


Fe / Ag(100) 


d Fe = 40 A 


0.66^ 




d Fc = 24 A 


0.65i& 




4 = 7 A 


2 3 46 




4 = 4 A 


5. 7 16 


Fe 4 / V 4 




1.25^ 


Fe 4 / V 2 




0.90^ 


Ni 


bulk 


2 ^34.47 


Cu / Co(lll) 




1.4&. 


Co / Cu(OOl) 


fee 


3.0^ 


Co 


fee, hard 


2.8±0.3 4 £ 




easy direction 


1.7 ±0.2^ 


Fe[001] 


bee 


0.0(63) 41 . 


Ni[001] 


fee 


0.(54)41 


Niflll] 


fee 


0.(45} 41 . 


Co [0001] 


hep 


0.0(36> 41 . 
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TABLE II: Third order Taylor series expansion coeffiecients of the total energy in case of 
Co/Cu n /Co. Notice that b 2 + 3ac > and (3c - a) 2 - 46 2 (not given here) are the smallest 
in magnitude for n = 30, see Eq. (|41l) . 



n 


a (meV) 


b (meV) 


c (meV) 


20 


0.18482 


1.21031 


1.34444 


21 


0.10742 


0.38908 


0.57940 


22 


-0.01466 


-0.46855 


-0.34401 


23 


0.05495 


0.33396 


0.38826 


24 


-0.01887 


-0.21263 


-0.19021 


25 


0.09246 


0.30579 


0.49206 


26 


0.06293 


0.65673 


0.56195 


27 


-0.05575 


-1.03111 


-0.78382 


28 


0.21518 


1.54457 


1.64362 


29 


-0.01959 


-0.25292 


-0.20745 


30 


-0.02395 


-0.17970 


-0.12111 


31 


0.20914 


1.64070 


1.60557 


32 


-0.24502 


-1.92659 


-1.93901 


33 


0.17413 


1.22408 


1.28963 
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TABLE III: Switching times (ns) as obtained by using the third order Taylor series expansion 
coemecients of the total energy, see Tab. [H] for Co/Cu n /Co m . FS (P or AP) denotes the final 
magnetic configuration by assuming that for the for the thick Co slab the direction of the moment 
is ni = +1. 



n 


m 


< 




FS 


7"min 


r(a bulk )3 Q 


r(a£ o/Cu/Co )£a 


r ( a scaled}44 


20 


13 


-1 


+1 


P 


0.01340 


1.33999 


0.54039 


0.108 24 


21 


12 


-1 


+1 


P 


0.02506 


2.50591 


1.01058 


0.186 97 


22 


14 


+1 


-1 


AP 


0.18977 


18.97722 


7.65309 


1.649 90 


23 


13 


-1 


+1 


P 


0.04415 


4.41567 


1.78074 


0.356 67 


24 


12 


+1 


-1 


AP 


0.13645 


13.64524 


5.50282 


1.01810 


25 


14 


-1 


+1 


P = AP 


0.02841 


2.84111 


1.14576 


0.24701 


26 


13 


-1 


+1 


P ~ GS 


0.05414 


5.41415 


2.18341 


0.43733 


27 


12 


+1 


-1 


AP 


0.06667 


6.66771 


2.68894 


0.49749 


28 


11 


-1 


+1 


P 


0.01154 


1.15445 


0.46556 


0.07902 


29 


13 


+1 


-1 


AP 


0.16889 


16.88991 


6.81132 


1.364 30 


30 


12 


-1 


+1 


P (AP) 


0.68836 


68.83831 


27.76096 


5.13610 


31 


14 


-1 


+1 


P 


0.01358 


1.35829 


0.54777 


0.11809 


32 


13 


+1 


-1 


AP 


0.01064 


1.06400 


0.42909 


0.08 594 


33 


12 


-1 


+1 


P 


0.01498 


1.49796 


0.60409 


0.11176 



23 



TABLE IV: Cross section of multilayer pillar sequence Co/Cu/Co used in experiments. 



multilayer pillar sequence 


cross section 


An (nm 2 ) 


Co(100 nm)/Cu(4 nm)/Co(d C o), 


5 4- 10 nm diameter 


(1^-4)- 19.63 


d Co = 2, 4, 7, 10 nmii 






Co(100 A)/Cu(60 A)/Co(25 A) 5 


130 i 30 nm diameter 


7853.98 4- 20106.19 


Co(40 nm)/Cu(6 nm)/Co(3 nm)£i 


~ 50 x 50 nm (sample 1) 


2500.00 




~ 130 x 60 nm (sample 2) 


7800.00 


Coftpiv^ nm.) /Cu(dr„ nm)/Coftpvon nm) 


< 100 nm diameter 


7853.98 


v r lxea — rree 






Co(30 nm)/Cu(10 nm)/Co(10 nm)S 


~ 60 -r 80 nm diameter 


2827.43 4- 5026.55 


Co(15 nm)/Cu(10 nm)/Co(2.5 nm)^ 


200 x 600 nm 2 


120000.00 


Co(10 nm)/Cu(10 nm)/Co(30 nm)S 


~ 40 nm diameter 


1256.64 


Co(3 nm)/Cu(10 nm)/Co(12 nm)ii 


~ 100 nm diameter 


7853.98 


Co(10 nm)/Cu(10 nm)/Co(30 nm^i 


~ 40 nm diameter 


1256.64 


Co(3 nm)/Cu(10 nm)/Co(12 nm)S 


0.05 x 0.10 /jm 2 


5000.00 




0.05 x 0.20 fj,m 2 


10000.00 




0.07 x 0.14 fim 2 


9800.00 




0.08 x 0.16 fim 2 


12800.00 


Co(3 nm)/Cu(10 nm)/Co(12 nm)^ 


50 -T- 200 nm circumference 


198.94 3183.10 


Co/Cu/Co(30 A)2 


0.05 x 0.10 fim 2 


5000.00 
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FIG. 1: C0/C1120/C0, spacer thickness: 34.65 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle 0. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle (top) and 
magnetoresistance as a function of the current (bottom), yj (A)) si / ( r min)si = 1> see Ea.(|44|). Solid 
lines serve as guidance for the eye. 
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FIG. 2: C0/C1121/C0, spacer thickness: 36.39 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle G (top) and 
magnetoresistance as a function of the current (bottom), (A)) si / ( r min)si = 1, see Eq . • Solid 
lines serve as guidance for the eye. 
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FIG. 3: C0/C1122/C0, spacer thickness: 38.12 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle G (top) and 
magnetoresistance as a function of the current (bottom), (A)) si / ( r min)si = 1, see Eq . • Solid 
lines serve as guidance for the eye. 
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FIG. 4: C0/C1123/C0, spacer thickness: 39.85 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle G (top) and 
magnetoresistance as a function of the current (bottom), (A)) si / ( r min)si = 1, see Eq . . Solid 
lines serve as guidance for the eye. 
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FIG. 5: C0/C1124/C0, spacer thickness: 41.58 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle G (top) and 



magnetoresistance as a function of the current (bottom), (^4o)si / ( r min)si = 1, see Eq, (f4"4"|) . Solid 
lines serve as guidance for the eye. 
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FIG. 6: C0/CU25/C0, spacer thickness: 43.32 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle G (top) and 
magnetoresistance as a function of the current (bottom), (A)) si / ( r min)si = 1, see Eq . . Solid 
lines serve as guidance for the eye. 
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FIG. 7: C0/CU26/C0, spacer thickness: 45.05 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle G (top) and 
magnetoresistance as a function of the current (bottom), (A)) si / ( r min)si = 1, see Eq . . Solid 
lines serve as guidance for the eye. 
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FIG. 8: C0/C1127/C0, spacer thickness: 46.78 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle Q (top) and 
magnetoresistance as a function of the current (bottom), (-Ao) SI / (T m i n ) SI = 1, see Ea.(|44|). Solid 
lines serve as guidance for the eye. 
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FIG. 9: C0/C1128/C0, spacer thickness: 48.51 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle Q (top) and 
magnetoresistance as a function of the current (bottom), y (A)) si / ( r min)gi = 1, see Ea.(|44|). Solid 
lines serve as guidance for the eye. 
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FIG. 10: C0/C1129/C0, spacer thickness: 50.24 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle G (top) and 
magnetoresistance as a function of the current (bottom), (A)) si / ( r min)si = 1, see Eq . . Solid 
lines serve as guidance for the eye. 
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FIG. 11: C0/C1130/C0, spacer thickness: 51.98 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle G (top) and 
magnetoresistance as a function of the current (bottom), (-Ao) SI / (T m i n ) SI = 1, see Eq . • Solid 
lines serve as guidance for the eye. 
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FIG. 12: C0/CU31/C0, spacer thickness: 53.71 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle G. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle Q (top) and 
magnetoresistance as a function of the current (bottom), yj (A)) si / ( r min)si = 1) see Ea.(|44|). Solid 
lines serve as guidance for the eye. 
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FIG. 13: C0/C1132/C0, spacer thickness: 55.45 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle 0. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle O (top) and 
magnetoresistance as a function of the current (bottom), */ (A)) si / ( r min)si = 1, see Ea.(|44|). Solid 
lines serve as guidance for the eye. 
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FIG. 14: C0/C1133/C0, spacer thickness: 57.18 A. Left column: twisting energy and sheet resistance 
as a function of the rotation angle 0. The dashed-dotted line refers to the first order approximation 
for the twisting energy. Right column: current as a function of the rotation angle (top) and 
magnetoresistance as a function of the current (bottom), (A)) si / ( r min)si = 1, see Eq . . Solid 
lines serve as guidance for the eye. 
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FIG. 15: Top: Twisting energy as a function of both rotation angles for the system with 25 spacer 
layers of Cu. Bottom: precessional energy at = 90°. 
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